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ABSTRACT 

Detailed information on the environmental con- 
ditions on Mars are very desirable for the design 
of photovoltaic systems for establishing outposts 
on the Martian surface. This paper addresses the 
variation of 'solar insolation (global, direct and 
diffuse) at the Viking lander's locations and can 
also be used, to a first approximation, for other 
latitudes. The radiation data is based on measured 
optical depth of the Martian atmosphere derived 
from images taken of the sun with a special diode 
on the Viking cameras; and computation based on 
multiple wavelength and multiple scattering of the 
solar radiation. The data are used to make esti- 
mates of photovoltaic system power, area and mass 
for a surface power system using regenerative fuel 
cells for storage and nighttime operation. 

INTRODUCTION 

NASA, through its Project Pathfinder, has put 
in place a wide-ranging set of advanced technology 
programs to address future needs of manned space 
exploration. Included in the mission sets under 
study is the establishment of outposts on the sur- 
face of Mars. The Surface Power Program in Path- 
finder is aimed at providing ultralightweight 
photovoltaic array technology for such an applica- 
tion (as well as for the lunar surface). Detailed 
information on the climatic conditions on Mars at 
the photovoltaic system location is very desirable. 
These include the variation of solar insolation; 
ambient temperature; albedo; and wind speeds and 
directions; with latitude, season, and time of the 
day. In addition, the effect of dust accumulation 
on the photovoltaic panels and radiation damage 
causing degradation of the output power needs to 
be assessed. 

As on the planet Earth, the solar radiation 
on the surface of Mars is composed of two compo- 
nents: the direct beam, and diffuse component. 

The direct beam is affected by scattering and 
absorption along the path from the top of the Mar- 
tian atmosphere to the Martian surface. Measure- 
ment of the optical depth [1,2] of the Martian 
atmosphere allows an estimate of the absorption 
and scattering out of the beam. These estimates 
were derived from images taken of the Sun and 
the Sun and Phobos with a special diode on the 
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Viking lander cameras. Since the cameras are sit- 
ting on the Martian surface, the measured intensity 
is directly related by Beer's law to the optical 
depth of the intervening atmospneric haze: 

G b = G 0 e * p [" i<5>] (,) 

where G 0 is the unattenuated irradiance at the 
top of the Martian atmosphere; Gb is the direct 
beam irradiance on Martian surface normal to the 
solar rays; x is the optical depth; and [m(z)] _1 
is the airmass determined by the zenith angle z. 

Earth-terresti al insolation data are accumu- 
lated over many years at different locations around 
the world and are given as long term average val- 
ues. The optical depth data for Mars are derived 
for less than two Mars years. Consequently, the 
calculated solar insolation corresponds to short 
term data. Furthermore , the measured opacities 
(optical depth) and the calculated irradiances per- 
tain to just two locations on the Mars planet; 
Viking lander 1 (VL1), located at 22.3° N 1 at i — 
tiude and 47.7° W longitude, and Viking lander 2 
(VL2) , located at 47.7° N latitude and 225.7° W 
longitude. However, the similarity in the proper- 
ties of the dust suspended above the two landing 
sites suggests that they are also representative 
of ones at other locations, at least, at latitudes 
not too far from the lander's sites. Data from 
lander VL1 may be used for latitudes 40° N to 40° S 
and data from lander VL2 for higher latitudes. The 
Martian atmosphere consists mainly of suspended 
dust particles, the amounts of which vary daily, 
seasonally, and annually, depending on local and 
global storm intensities and their duration. The 
optical depth values given in the section entitled 
OPTICAL DEPTH are assumed to be constant throughout 
the day. Large values of optical depth correspond 
to global storms, i.e., days with low insolation 
("dark days"). 

The albedo of the Martian surface varies in 
the range of about 0.1 to 0.4. The irradiances 
derived in the section entitled SOLAR RADIATION 
correspond to 0.1 albedo but can be also used for 
other values of albedo, to the first 
approximation. 

In this paper we calculated the variation of 
solar irradiance which is the major environmental 
component for photovoltaic system design. We do 
not report here data on ambient temperature, wind 
speed and direction, albedo and the effect of dust 
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New information about Mars climate may be forth- 
coming in the future based on new analysis, obser- 
vations and future flight missions. The Mars 
climatic data for photovoltaic system design will 
thus be updated accordingly. 

We have also completed a preliminary calcula- 
tion of photovoltaic power system size and mass 
(including storage), based on the data now in hand, 
for use on the Martian surface during a so-called 
"sprint" mission (i.e., 30 day duration or less on 
the surface). The detailed assumptions appear in 
the last section of the paper on power system con- 
siderations, and the results appear in Table II. 

OPTICAL DEPTH 

The most direct and probably reliable esti- 
mates of opacity are those derived from Viking 
lander imaging of the Sun. Figures 1 and 2 show 
the seasonal variation of the normal -i nc i dence of 
the optical depth at the Viking lander locations 
VL1 and VL2 , respectively. The season is indi- 
cated by the value of L s , areocentric longitude 
of the Sun, measured in the orbital plane of the 
planet from its vernal equinox (L s = 0° and 180° 
corresponding to spring and fall equinox for the 
northern hemisphere, respectively; and L s = 90 
and 270° corresponding to northern and southern 
summer solstices, respectively). Figures 1 and 2 
were derived from references by Pollack [1,2] and 
Zurek [3] and were discretized for each 5°. As 
mentioned before, the optical depth is assumed to 
remain constant throughout the day. Opacities are 
minimum during the northern spring <L S = 0° to 90°) 
and summer (L s = 90° to 180°), and maximum during 
southern spring (L s = 180° to 270°) and summer 
(L s = 270° to 360°), the seasons during which most 
local and major dust storms occur. When dust 
storms are not present, the optical depth are typi- 
cally about 0.5. Two global dust storms occurred 
during the period of the observation as indicated 
by the high values of the optical depth (lower 
bound values) . 

GLOBAL AND LOCAL DUST STORMS 

The intensity of Martian global and local dust 
storms is defined in terms of opacity of the dust 
it raises. Global dust storms are those which 
obscure planetary-scale sections of the Martian 
surface for many Martian days (sols), whereas local 
dust storms are less intense, and form and dissi- 
pate in a few days or less. From photovoltaic sys- 
tem design point of view, the intensity, frequency, 
and duration of these storms may be viewed as "par- 
tially cloudy" and "cloudy" days for which addi- 
tional energy storage in the photovoltaic system 
must be taken into account. The characteristics 
of global and local dust storms are listed below. 

Global Dust Storms 

(1) One, or occasionally two global dust 
storms of planetary scale may occur each Martian 
year. The duration may vary from 35 to 70 days or 
more. Although global dust storms do not occur 
every year, their occurrence is fairly frequent. 


^2) Global dust storms begin near oerihelion, 
when solar insolation is maximum 'southern spring 
and summer) in the southern mid-lari tude . 

'3) The first global dust storm observed by 
VL (1977) spread from a latitude of 40° S to a lat- 
itude 48° N in about 5 to 6 days. 

(4) The opacity during the global dust storm 
i s greater than 1 . 

Local Dust Storms 

(1) Local dust storms occur at almost all lat- 
itudes and throughout the year. However, they 
have been observed to occur most frequently in the 
approximate latitude belt 10° to 20° N and 20° to 
40° S, with more dust clouds seen in the south 
than in the north, the majority of which occurred 
during southern spring. 

(2) Based on Viking orbiter observations, it 
is estimated that approximately 100 local storms 
occur in a given Martian year. 

(3) Local dust storms last a few days. 

(4) The opacity of local dust storms may be 
assumed about 1 . 

SOLAR RADIATION 

"he solar insolation on the surface of Mars 
is composed of the direct beam and the diffuse com- 
ponents. The net solar flux integrated over the 
solar spectrum, on the Martian surface was calcula- 
ted by Pollack [53 based on multiple wavelength 
and multiple scattering of the solar radiation. 
Derived data of this calculation is shown in 
Table I by the "normalized net flux function" 
f(z,x) where the parameters are the zenith angle 
z and the ootical depth t. This table pertains 
to albedo of 0.1 but can be used for higher albedo 
values to a first approximation. Using these data 
we calculated the global solar irradiance. We 
assumed that the diffuse irradiance is obtained by 
subtracting the beam from the global irradiance. 

The solar irradiance components, on a horizon- 
tal Martian surface, are related by: 

G h = Gb h + G<j|, (2) 

where 

Gh global irradiance on a horizontal surface 
Gbh direct beam irradiance on a horizontal surface 
G<jh diffuse irradiance on a horizontal surface 

The diffuse irradiance of the Martian atmospnere 
may be a result of a different mechanism than that 
for the Earth atmosphere, nevertheless, to a first 
approximation, we will apply Eq. (2) as for Earth- 
terrestrial calculations, the beam irradiance 
Gbh on a horizontal surface is given by: 
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For zenith angles z < 70° one can approximate the 
airmass by: 

1 a _] (4) 

m(z) * cos z 


The beam irradiance at the top of the Martian 
atmosphere is given by: 



(5) 


where S i s the solar constant at the mean Sun- 
Earth distance of I AU, i.e., S * 1371 W/m 2 ; r is 
the instantaneous Sun-Mars distance (heliocentric 
distance) given by C6J: 


_ a(l - e 2 ) 

1 + e cos 9 


( 6 ) 


where a is the Mars semimajor axis, and e is 
the Mars eccentricity, i.e., e * 0.093377; and 9 
is the true anomaly given by: 

9 = L $ - 248° (7) 

where L s is the areocentric longitude and 248° is 
the areocentric longitude of Mars perihelion. The 
Sun-Mars mean distance in AU units is 1.5236915, 
therefore, the mean beam irradiance at the top of 
Mars atmosphere is: 1371/1 . 523691 5 2 = 590 W/m 2 . 

The instantaneous beam irradiance is given by 
Eqs . (5) to (7): 


590 


^1 + e cos(L - 248° ) j 


( 3 ) 


where w is measured from the true noon westward. 

The following figures of solar irradiances 
were calculated based on Table I data and the mean 

irradiance of 590 W/m 2 . Figures 6 to 8 describe 
the variation of the global, beam and diffuse irra- 
diances, respectively, on a horizontal Martian sur- 
face; and are given in pairs as function of the 
optical depth t, and zenith angle z. The varia- 
tion of the global irradiance , Eq. (9), is 
shown in Figs. 6(a) and <b). The beam irradiance 
Gph is obtained using Eq. (3) and is shown in 
Figs. 7(a) and (b). The beam irradiance shows a 
sharp decrease with increasing optical depth, and 
a relative moderate decrease with increasing 
zenith angle. The diffuse irradiance G<jh is 
obtained from Eqs. (2), (3) and (9) and is shown 
in Figs. 8(a) and (b). The diffuse irradiance 
shows a sliding maximum with the variation of the 
zenith angle. 

For a given day l s and latitude <$>, one can 
calculate the variation of the zenith angle z as 
function of the solar time T using Eqs. (10) to 
(12). Referring to Figs. 1 and 2 for the given 
l s , the optical depth x is determined, and using 
Figs. 6 to 8 one can calculate the variation of the 
solar irradiance for the given day. The results of 
such a procedure is shown, for example, in Fig. 9 
where the variation of the global, beam and diffuse 
irradiance on a horizontal surface is depicted as 
function of the solar time for L s * 153° and at 
Viking lander VL1 . 

POWER SYSTEM CONSIDERATIONS FOR MARTIAN 
SURFACE OPERATION 


and is shown in Fig. 3. 

The global irradiance Gq on a horizontal 
Martian surface is obtained by: 


The cosine of the zenith angle is: 

cos z = sin 4> sin 6 + cos cos 8 cos w (10) 
where 

$ latitude 
8 dec 1 i nation angle 
a) hour angle 

The solar declination angle is given by: 

sin 8 * sin 8 q sin (11) 

where 8 0 = 24.936° is the Mars obliquity of rota- 
tion axis. The variation of the solar declination 
angle is shown in Fig. 4. 

The relation between the planet solar time T and 
the hour angle u, is shown in Fig. 5 and is given 
by: 

CO =* 1 51 - 180 (12) 


The data presented in this paper represent a 
start on a longer term effort to provide a more 
complete analysis of temperature and insolation 
levels on Mars. The purpose for doing so is two- 
fold: (1) to establish the feasibility of using 

photovoltaic power systems on the surface of Mars, 
and (2) to allow for much more accurate calcula- 
tions of critical power system character! st i cs , 
such as area and mass, than has been possible in 
the past. Of major concern are the dust storms, 
which have been observed to occur on local as well 
as global scales, and their effect on solar array 
output. In general, the assumption has been that 
global storms would reduce solar array output 
essentially to zero, because the opacity of the 
atmosphere may rise to values ranging from 3 to 9, 
depending on the severity of the storm. Further- 
more, such high values of opacity may persist for 
long periods of time (often several weeks), so 
that predicted requirements for energy storage 
quickly become much too large to be practical. As 
shown in the preceding discussion, however, the 
opacity primarily affects the direct beam compo- 
nent of the solar radiation. Hence there is still 
an appreciably large diffuse component, even at an 
atmospheric opacity of 4 (Fig. 8), with the result 
that solar array operation is still possible. 
According to the data shown in the figure, array 
output at midday could be expected to be over 
40 percent of its “clear weather'* midday output 
during the height of that particular storm. 
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The second, and primary use for data such as 
are presented here, is to allow for more accurate 
estimates of photovoltaic power system sizes and 
masses in system analysis and trade-off studies. 
Ultimately, of course, such data will be required 
for final engineering design should a manned mis- 
sion to Mars actually be attempted at some future 
date. Prior to any such actual engineering phase, 
however, it is important to analyze and compare 
all the relevant technology options, and to select 
those which look most favorable for further study 
and development. As an example, a comparison will 
be made of array areas and masses for two rad- 
ically different solar cell options. It is 
assumed that energy storage will be provided by 
regenerative f oe 1 cells (RFC) using hydrogen and 
oxygen as reactants. The storage system efficien- 
cy in this calculation takes into account the 
storage "round- tr i p" efficiency as well as the 
power management and distribution efficiencies. 

For storage times on the order of 10 hr or more, 
RFC energy storage densities have been shown to 
approach 1000 W-hr/kg [7]. The two solar cell 
options are thin film amorphous silicon at 
12 percent AMO efficiency, and thin, single crys- 
tal GaAs cells at 26 percent AMO. Both cell 
types are assumed to be integrated into a suitable 
iightweight, flexible blanket assembly. Based on 
projections from the NASA Advanced Photovoltaic 
Solar Array Study [8], the GaAs blanket specific 
mass is calculated to be 0.5 kg/m 2 and the aSi 
olanket specific mass is 0.04 kg/m 2 . The array 
structure in both cases is assumed to have a mass 
of 0.56 kg/m 2 meter of array area, to which must 
oe added the blanket mass for each case. No tem- 
perature effects are included in this calculation, 
and it is assumed that the efficiencies quoted are 
at the expected operating temperature of the 
array. Also, not included is the mass of the 
power conditioning and distribution equipment for 
the loads on the power system. 

The insolation level and opacity data are 
that observed by VL1 on day L s = 153V The actual 
irradiance level of 590 W/m 2 during this portion 
of Mar's orbit is approximately equal to the orbi- 
tal mean value, and remains within 10 percent of 
that value for several weeks on either side of 
that day. It is also a period of relatively sta- 
ble atmospheric conditions, with opacities around 
0.5 for several days at a time. Figure 9 shows 
the variation with zenith angle of the irradiance 
on a fixed horizontal surface for that same day. 
Table II compares the area and mass of the two 
solar array types when both are fixed and 
horizontal, as might be the case if operational 
simplicity is desired for their deployment. 

A rray power per unit area and unit mass have 
been calculated with the condition that the array 
must provide its full rated output of 25 kWe for 
at least 5.5 hr past Martian noon, at which time 
the solar irradiance is 70 W/m 2 on the day selec- 
ted. The ratio of the integrated total energy to 
the peaK solar insolation is 0.59, which yields a 
daily average irradiance of 0.348 kW/m^. It is 
interesting to note that although the mass of the 
aSi blanket was calculated to be less than 
10 percent of the GaAs blanket, the aSi system 


mass is nearly 20 percent greater than :ne GaAs 
system mass. The reason is, of course, that me 
lower efficiency projected for otSi results in a 
much larger total area, and hence total mass, Vr 
that array, under the assumption that the same 
structural mass is required in each case. Clearly, 
a more detailed design is required to resolve the 
issue. 

CONCLUSIONS 

The major environmental parameter needed for 
photovoltaic system design is the variation of the 
solar insolation. The data and procedure from 
which the solar insolation can be determined is 
given in this paper, and is based on the normalized 
net solar flux function f(z,x>. The opacities 
were measured at the two Viking lander locations 
but can also be used, to the first approx i ma t i on . 
for other locations. Figures 6 to 9 show the glo- 
bal, beam and diffuse irradiances on the Horizon- 
tal surface. Any desired insolation quantities on 
any inclined surface can be further derived based 
on the f ( z , t ) function and the appropriate sc’ar 
calculations. 

An example is presented for a surface power 
system using regenerative fuel cells for storage, 
and two different array types--one with GaAs cells, 
and one with aSi thin film cells. The results 
demonstrate the ease with which the solar 
insolation data provided here can be used to 
determine inportant power system parameters for 
system comparison studies, and for final 
engineering design. 
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TABLE I. - NORMALIZED NET FLUX FUNCTION f(2,t) AT THE MARTIAN SURFACE 



TABLE II. - MARS MISSION POWER SYSTEM COMPARISON 

Total power delivered to load (kWe): 25 

Solar irrandance (kW/m 2 ): 0.348 (daily ave.). 


(a) Performance parameters 


Technology 


Array 


RFC energy 


System total 

Array 





storage mass, 
*9 





sized 

needed, 

KW 


•^/kHe 

kg/kHe 

Mass-kg 

Area-m 2 

Mass-kg 

kg/kHe 

Mars/Ground/GaAs/RFC 

65.0 

69.1 

1728 

516 

1624 


2244 

H9 

■m 

Mars/Ground/aSl/RFC 

140.7 

84.3 

2107 

527 

3518 


2634 

Baa 

Kali 


(b) Input parameters 
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FIGURE 1. - OPTICAL DEPTH AS MEASURED FOR VIKING LANDER VL1 
AS FUNCTION OF AREOCENTRIC LONGITUDE. 




















SOLAR IRRADIANCE AT TOP OF MARS ATMOSPHERE, OPTICAL DEPTH 
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FIGURE 2. - OPTICAL DEPTH AS MEASURED FOR VIKING LANDER VL2 
AS FUNCTION OF AREOCENTRIC LONGITUDE. 
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FIGURE 3. - SOLAR IRRADIANCE AT TOP OF MARS ATMOSPHERE AS 
FUNCTION OF AREOCENTRIC LONGITUDE. 
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FIGURE *1. - VARIATION OF SOLAR DECLINATION 
ANGLE 6, WITH AREOCENTRIC LONGITUDE, L $ . 
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FIGURE 5. - SOLAR TIME AND HOUR ANGLE RELATION. 
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GLOBAL IRRADIANCE, G h , W/m 
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OPTICAL DEPTH, T 

(A) EFFECT OF OPTICAL DEPTH WITH SUN ZENTH ANGLE AS A 
PARAMETER. 



(B) EFFECT OF SUN ZENITH ANGLE WITH OPTICAL DEPTH AS A 
PARAMETER. 

FIGURE 6. - VARIATION OF GLOBAL IRRADIANCE WITH OPTICAL DEPTH 
AND SUN ZENITH ANGLE ON A HORIZONTAL SURFACE. 
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(A) EFFECT OF OPTICAL DEPTH WITH SUN ZENITH AN6LE AS A 
PARAMETER. 
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FIGURE 7. - VARIATION OF BEAM IRRADIANCE WITH OPTICAL DEPTH 
AND SUN ZENITH ANGLE ON A HORIZONTAL SURFACE. 
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(B) EFFECT OF SUN ZENITH ANGLE WITH OPTICAL DEPTH AS A 
PARAMETER. 


9 





DIFFUSE IRRADIANCE, G hh , W/m' 
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OPTICAL DEPTH/ T 

(A) EFFECT OF OPTICAL DEPTH WITH SUN ZENITH ANGLE AS A 
PARAMETER. 
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(B) EFFECT OF SUN ZENITH ANGLE WITH OPTICAL DEPTH AS A 
PARAMETER. 

FIGURE 8. - VARIATION OF DIFFUSE IRRADIANCE WITH OPTICAL 
DEPTH AND SUN ZENITH ANGLE ON A HORIZONTAL SURFACE. 
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SOLAR TIME 

FIGURE 9. - DIURNAL VARIATION OF SOLAR IRRADIANCE ON A 
HORIZONTAL SURFACE AT VIKING LANDER VL1 . AREOCENTRIC 
LONGITUDAL, 153; OPTICAL DEPTH, 0.5. 
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